Point-contact tunneling involving low-dimensional spin-triplet superconductors 
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We modify and extend previous microscopic calculations of tunneling in superconducting junctions 
based on a non-equilibrium Green function formalism to include the case of spin-triplet pairing. We 
show that distinctive features are present in the I-V characteristics of different kinds of junctions, 
in particular when the effects of magnetic fields are taken into account, that permit to identify the 
type of pairing. We discuss the relevance of these results in the context of quasi one-dimensional 
organic superconductors like (TMTSF)2PF6 and layered compounds like S^RuCU. 
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Unconventional superconductivity is a generic classi- 
fication for all those superconducting materials that do 
not fit into the s-wave spin-singlet pairing picture of stan- 
dard BCS theory[l|. One of the possible extensions of 
that theory considers anisotropic pairing mechanisms ac- 
companied by spin-triplet states. One such unconven- 
tional scenario is realized in the p-wave spin-triplet su- 
perfiuid state of 3 He. Among the superconductors the 
appearance of unconventional superconductivity, always 
accompanied by strong correlations, results in most cases 
in spin-singlet pairing. The search for triplet pairing su- 
perconductivity in strongly correlated electron systems 
becomes thus natural. It was proposed^, for instance, 
in the case of Sr2Ru04 and has to date a consider- 
able amount of experimental support This compound 
shows, in common with the high- T c cuprates, a layered 
structure. Other candidates for triplet pairing, though 
the evidence is less conclusive, are in the class of so called 
organic superconductors 3 like f.i. (TMTSF) 2 X, with X 
a counter anion as PFg, CIO4, etc. Materials of this class 
are low-dimensional conductors in their normal state. 

A distinctive feature of triplet superconductors is the 
possibility of large upper critical fields; this is contrary 
to BCS superconductors for which the paramagnetic re- 
sponse of the normal state puts a limit to the field up to 
which the superconducting phase is energetically favor- 
able even in the absence of orbital effects (the Clogston 
or BCS Pauli paramagnetic limit, H p ). It was indeed ob- 
served during the last few years that at low temperatures 
the above mentioned organic compounds show upper crit- 
ical fields that exceed H p , usually depending on the an- 
gular orientation of the field with respect to the crys- 
talline axesQ- In general, for layered and chained com- 
pounds, superconductivity at high fields is attained when 
the field is oriented so that orbital effects are suppressed. 
Although high critical fields constitute an indication of 
triplet-pairing, unambiguous interpretation of this effect 
is difficult because one needs to rule out other possible 
phases able to accommodate large magnetic fields (e.g. 
the Larkin-Ovchinikov-Fulde-Ferrel phase) 0. 

It would be highly desirable to have a direct indication 
of the nature of the pairing, as in the case of cuprates. A 
possible experimental determination of the symmetry of 
the pairing could come from tunneling experiments. In 



the past some of the most crucial experimental verifica- 
tions of BCS theory came from tunneling experiments; 
and, after the discovery of the Josephson effect, some 
of the most important practical applications of super- 
conductivity involve tunneling junctions. In addition, 
tunneling per se acquired renewed relevance with the 
development of scanning tunneling microscopy (STM), 
actively used nowadays in the study of superconductiv- 
ity. The simplest theoretical models used to interpret 
tunneling experiments in superconductors are those that 
go under the name of semiconducting band models^, 0|. 
A more systematic approach is that based on tunnel- 
ing Hamiltoniansjil lit fiof. Two limiting cases are usu- 
ally considered in the calculations: planar interfaces and 
point contacts (the latter usually realized in experiments 
using break junctions, pressed crossed wires or close STM 
contacts) . 

In this letter we focus our attention on the point- 
contact case. We extend previous tunneling Hamiltonian 
calculations to include the case of triplet-pairing super- 
conductors. As compared to previous calculations, we 
simplify considerably the formalism making it more ver- 
satile and easy to implement. Our aim is to explore the 
possible use of tunneling experiments as a way of charac- 
terizing the type of pairing in unconventional supercon- 
ductors. Using non-equilibrium Keldysh Green functions 
we calculate the full current- voltage characteristics of dif- 
ferent types of tunnel junctions mixing normal metals, 
singlet and triplet superconductors. We also study the 
effects of applied external magnetic fields on the trans- 
port properties of the different junctions. 

Our starting point is the tunneling Hamiltonian, 



-ffjun — Hi + H 2 + H tun 

H tnn = -t £ CT U\ (x = 0) ^ _ (x = 0) + h.c. 



(f) 



The first two terms describe the two leads (supercon- 
ducting or otherwise) and the third one models the tun- 
neling processes in which an electron with spin a hops 
from one lead into the other. The current is propor- 
tional to the rate of change in the relative particle num- 
ber, I = f (d t (N 2 - Nx)) = ^ ([H tun , Ni - A/2]). 

To model the leads in calculations of point-contact 
transport on conventional superconductors, simple mod- 
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els suffice to achieve quantitative agreement with the 
experiment. Since dimensionality plays no role in the 
tunneling, contrary to the case in some planar junction 
experiments, all the standard calculations can be formu- 
lated using one-dimensional leads. The situation is more 
complex for unconventional superconductors, for which, 
the anisotropic nature of the pairing has to be taken into 
account. However, the ruthenates and organic supercon- 
ductors are deemed to have p-wave symmetry, and since 
both s-wave and p-wave symmetries can be modeled in 
a one dimensional chain, we can set up a formalism that 
comprises both cases. We consider a one dimensional 
band with two Fermi points and expand the fermion fields 
around them defining left (L) and right (R) moving fields. 
Using these fields, and in the spirit of BCS theory, we in- 
troduce the four gap functions: 

A a (x) oc (a iP L - (x) G a a(i (a)) (2) 

where a, (3 £ (|,T) = ( — 1 , + 1 ) are summed over (a = 
—a); a = 0, . . . , 3; <r°o is the identity matrix while the 
other three are the usual Pauli matrices. With this def- 
inition Ao (x) is the conventional spin-singlet order pa- 
rameter and the other three gap functions form a vector 
of spin-triplet order parameters |3J, A(x) — A(x)d(x). 
Approximating the order parameter to have no spatial 
dependence, we write in Fourier space 



K n — Qka^cka^cka 



A 



Rkf3 u /3a 



a 1p 



Lki 



h.C 



where K n = H n — ^, n N n with \i n the chemical potential 
of lead n. All the indexes are summed over, in particular 
c G {L,R) = (— 1,+1) sums over the two possible chi- 
ralitics and £™ fcQ = ck — fi n — ah are the corresponding 
linear dispersions, shifted by the inclusion of chemical 
potential and magnetic field along the z-axis. We take 
the quantization axis ( z ) along the field direction and 
consider the cases of triplet order parameters parallel or 
perpendicular to it. 

Since we are dealing with an out of equilibrium phe- 
nomenon, we use Keldysh formalism|llj to treat the tun- 
neling term to all orders, calculating the full I-V line 
and giving a quantitative account of its subgap struc- 
ture. One past implementation of such an approach, for 
the s-wave case, reduced the problem to the solution of 
a set of linear recursion relations Here, instead, we 
notice that in this formalism the current becomes 



7 = | E/ -^{4,^-4^ 2 ,c 



kel 



(3) 



where 'kel' denotes the Keldysh component of the corre- 
lation function and the ij) stand for ^(uj,x = 0). Since 
the current depends only on the fields at x — one can 
analytically integrate the x dependence in the leads to 
obtain from Eq. a local and quadratic Keldysh ac- 



Here Stun is obtained from fltun and S n is the local ac- 
tion of lead n, of the form S n = J 
Where ty n is an 8 component spinor (a Keldysh extended 
Nambu-Eliashberg spinor), and g~ x is a matrix whose 
components can be computed from K n . Its inverse (g) is 
given by the standard advanced, retarded and Keldysh 
components of the local Green functions of the lead. As 
an example we give the expressions for the case when 
Ai = A 2 = 0, 



[ret , ad v] 
ilc<7,ca 



[ret, ad v] 
9ca,c<7 



(u> — fj, n + coh ± i0 + ) 



2 J |A + co- A 3 | 2 - (w - /Lt„ + coh ± iO+f 



(A + caA 3 ) 1 



2 J |A + co- A 3 | 2 - (w - /Lt„ + coh ± iO+f 



The notation in the numerator of the anomalous func- 
tions means that complex conjugation is in order when 
c = L. The Keldysh component is obtained immediately 
as 5 kcl = ( 5 rct - ff adv ) tanh ((w - fx n ) /2T). 

Since the action is quadratic, it can be diagonalized 
which allows to compute the current from Jjyl. For nor- 
mal leads this can be done analytically. However, for 
superconducting leads, special attention must be paid to 
the fact that frequencies have different reference Fermi 
levels in each lead when there is a bias applied. Within 
each lead, frequencies with equal positive and negative 
shifts from the Fermi level form part of the same spinor 
element and the corresponding states are related by the 
coherent pairing processes in the superconductors; across 
leads, same frequency states are related by the tunneling 
matrix elements of the action. As a result, the full action 
for the junction is no longer frequency diagonal. How- 
ever, since it is quadratic it can still be written explicitly 
as a matrix. To each value (o>o) m the frequency window 
(of size eV = [ii — /12) defined by the chemical potentials 
in the two leads, one infinite set of related frequencies 
can thus be assigned (p > 0): 



2/i2-pmod2 — ^p-l 
2/il+pmod2 — ^l-p 



(4) 



tion for the fields at x = 0: 



Si 



These sets are all independent and the action is block di- 
agonal between different ones. Discrctizing the frequen- 
cies in this window, automatically defines a discretization 
of the whole frequency axis. This allows for a numerical 
solution of the problem: we deal with one set of frequen- 
cies at a time, and since the sets are infinite, we trun- 
cate their hierarchy at some distance from the central 
frequency window. This is equivalent to introducing a 
soft limit in the number of allowed Andreev reflections 
(Aa). We can then numerically invert the corresponding 
block-diagonal action, written as a matrix in frequency 
space. The off-diagonal Green function matrix elements 
thus obtained allow to compute the current using Eq. 

The practical implementation of this approach is quite 
simple and allows to consider the (combined) effects of 
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finite temperature, applied magnetic fields, contact po- 
tentials in the junction, spin-flip tunneling or spin-flip 
scattering processes in the leads. It is also possible to 
compute the a.c. response. Our interest is in compar- 
ing singlet and triplet superconductor junctions and how 
they respond differently in the presence of an external 
field; we will ignore other additional complications. 

In order to illustrate the different I-V characteristics 
for different types of junctions, we fix a set of parameters. 
We take t = 0.2 for the tunneling overlap integral, and 
when there is a magnetic field we fix its value to h = 0.2 
in units of A (the magnitud of the singlet gap, A , or 
of the triplet vector order parameter depending on the 
case). All the curves we show are for the d.c. response in 
the limit of vanishing temperatures. For the truncation 
procedure we take — 3 (and verify that larger values 
produce identical curves). Let us use the notation, N: 
normal-metal, S: singlet-superconductor and T: triplet- 
superconductor. We show in Fig.Q-(a) typical curves for 
an N-S junction. The effect of the field (for any orien- 
tation) is to produce what would be seen as a Zeeman 
splitting in the differential conductance peak. Notice the 
sub-gap shoulder on the I-V curve when eV < A (for 
h = 0); its origin is in the coherent Andreev processes 
that take place in the junction contact. 

We show in Fig.EJ(b) typical curves for an N-T junc- 
tion. The solid line corresponds to zero field and the 
dashed line is for that same junction but in the presence 
of a field that is aligned with the vector order parame- 
ter. If one considers a field perpendicular to the order 
parameter (h _L d), the I-V characteristic remains unaf- 
fected, i.e. identical to the one for zero field. Notice the 
absence of a sub-gap shoulder on the I-V curve, a result 
of the odd real-space symmetry of the superconductor 
probed locally by an ideal point-contact junction. Given 
the model we consider, these results would be relevant for 
tunneling experiments into a sample edge parallel to the 
chains, when a zero-bias peak is not expected ^lEH- 

Let us now examine the case of S-T junctions, see 
Fig. EJ(c). One dotted line is the N-N characteristic 
as in the other cases; the other is the I-V curve of an 
S-S junction that shows all the features well known in 
the literature[3 E3, in particular a 'sub-gap' shoulder 
with Andreev steps at eV = 2A/n (with n = 1,2,...). 
This curve is, bar orbital effects, not sensitive to applied 
fields. The solid line is always the same regardless of the 
orientation of the vector order parameter on the triplet- 
pairing side and the current amplitude is systematically 
smaller than in the S-S case. Distinctively, the 'sub-gap' 
structure shows only the first two steps and the current 
becomes zero for eV < A. Regarding the effect of an 
applied field, the curve remains unchanged if the field 
is parallel to the vector-order-parameter direction, but 
shows a Zeeman effect if the field is perpendicular to it 
(dashed line) . This is to be confronted with the N-T case 
where the splitting accompanies a field h || d. 
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FIG. 1: Three different sets of I-V characteristics vertically 
displaced for clarity. From top to bottom: (a) for an N- 
S junction, (b) N-T junction, (c) S-T junction. The solid 
(dashed) lines correspond to the characteristics in the absence 
(presence) of an oriented magnetic field (see text for details) . 
In all three sets is plotted as reference the curve for the N-N 
case (diagonal straight dotted lines) and in the third one the 
curve for an S-S junction is also given. 



Let us make contact with the experimental situation on 
the different compounds afore mentioned, and dwell first 
on the particularly interesting case of (TMTSF^PFg. 
This quasi-one-dimensional organic compound shows the 
largest deviations from the H p limit among the different 
members of the Bechgaard salts family, and also as com- 
pared with other organic superconducting salts A re- 
cent series of experiments focused on measuring H C 2 (T) 
for this material[H H H H E3 . The observed large 
deviations from the H p limit as well as the anisotropy 
and angular dependence of the upper critical field make 
the case for an equal-spin triplet-pairing superconduct- 
ing phase with an order parameter oriented mainly along 
the chains (lij. Other evidences for spin-triplet pairing 
in (TMTSF) 2 PF 6 also exist, but as with the H c2 mea- 
surements, most of these experiments do not probe the 
spin parity directly and their conclusions are difficult to 
interpret due to the uncertainties regarding the orbital 
symmetry of the superconducting phase; however, re- 
cent NMR Knight Shift measurements add strong sup- 
port to the spin-triplet scenario [l9j. Here we argue that 
the magnetic field response in N-T point-contact tunnel- 
ing experiments would also constitute a direct probe of 
the spin-pairing symmetry (in contrast with tunneling 
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across planar junctions that is sensitive to the orbital- 
pairing symmetry[20j). H C 2 measurements show that 
for magnetic fields along the direction of the conduct- 
ing chains (a-axis) the upper critical field is paramag- 
netically limited|l8|]. This corresponds, in our conven- 
tion, to a vector order parameter aligned with the field 
(h || d). In this case a Zeeman splitting of the differential 
conductance peak, similar to that in conventional super- 
conductors (N-S junctions), should be observed. As the 
field is rotated, the splitting would be suppressed and 
for a magnetic field oriented parallel to the b'-axis there 
would be no Zeeman effect (accompanied by the possibil- 
ity of applying large fields that are not paramagnetically 
limited). The disappearance of splitting even as the field 
is being increased would constitute a clear signature of 
spin-triplet superconductivity. 

As in the case of N-T junctions, we can envisage using 
the Zeeman response of S-T junctions as a direct probe 
for spin-triplet order. If for instance, a field is applied 
along the b'-axis of (TMTSF) 2 PF 6 , we predict a Zee- 
man splitting of the main differential conductance peak. 
This would constitute a clear sign of unconventional su- 
perconductivity since such an effect does not take place 
for standard BCS superconductors (S-S junctions). The 
b' direction is the one on which the upper critical field 
is not paramagnetically limited, so relatively large fields 
could be applied in order to obtain a clear signal. 

Our considerations could be extended to the case 
of the layered compounds believed to be triplet 
superconductors 0. Among them, Sr2Ru04 is the best 
studied so far , bu t only few tunneling experiments were 
performed[2ll HlEl I2I EoT] . and none so far with good 
resolution in the presence of an applied external mag- 
netic field to observe Zeeman effects. One of the con- 
spicuous features observed in some of these experiments 
is the presence of a zero-bias anomaly in the differen- 
tial conductance. Its explanation is still a matter of de- 
bate, but seems to require extended contacts and sign 
changing order parameters. To include the effect of 
'zero energy states' at the interfaces would require ex- 
tensions to our scheme, possibly incorporating certain 
aspects of those calculations already done for planar 
i unctions [I21 I2H 121I Ejj . Our general findings about the 
effect of magnetic fields should however apply, since they 
refer to features to be measured at voltages of the order 
of the superconducting gap. 

Summarizing, the point-contact tunneling involving 
unconventional superconductors with spin-triplet pairing 
displays interesting characteristic features. In particu- 
lar the Zeeman response to an external magnetic field 
is such that it may allow for the identification of triplet 
phases and might be relevant for future experiments. The 
prediction of a truncated sub-gap structure in S-T junc- 
tions is also very interesting. For future extensions of 
this work, different effects (finite temperature, contact 
potentials in the junction, spin-flip tunneling, etc.) can 



be easily incorporated to our calculations. 
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